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Recent results from the Apollo Seismic Network suggest that (1) primitive differentiation 
occurred in the outer shell of the Moon to a depth of approximately 300 km and (2) the 
central region of the Moon is presently molten to a radius of between 200 and 300 km. If 
early melting to a depth of 300 to 400 km was a consequence of accretional energy, very 
short accretion times are required. The best model for the zone of original differentiation 
appears to be a crust 40 to 80 km thick, ranging in composition from anorthositic gabbro 
to gabbro, and overlying an ultramafic cumulate (olivine-pyroxene) about 250 km thick. 
The best candidate for the molten core appears to be iron or iron sulphide. A new class of 
seismic signals has recently been identified that may correspond to shallow moonquakes. 
These are rare, but much more energetic than the more numerous, deep moonquakes. 


Five seismic stations have been placed on 
the Moon by the astronauts of Apollo mis- 
sions 11, 12, 14, 15, and 16. The Apollo 11 
station, powered by solar cells and intended 
for operation only during the lunar day, 
failed after exposure to the first nighttime 
period. The remaining four stations, pow- 
ered by radio isotope thermal generators 
(RTG), have operated continuously since 
their initial activation. These four stations 
constitute the Apollo Seismic Network. The 
locations and installation dates of these sta- 
tions are given in table 1. 

The network was completed as of April 
1974, with the installation of the fourth sta- 
tion in the Descartes region of the southern 
highlands during mission 16. It is estimated 
that the maximum lifetime of each station 
will be about 10 yr, corresponding to the use- 
ful lifetime of the RTG power source. Each 
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station contains four seismometers. Three of 
these seismometers form a triaxial set (one 
sensitive to vertical motion and two sensitive 
to horizontal motion), with sensitivity to 
ground motion sharply peaked at .45 Hz. The 
fourth seismometer is sensitive to vertical 
motion, with peak sensitivity at 8 Hz. These 
instruments can detect vibrations of the lu- 
nar surface as small as 1/^ A at maximum 
sensitivity. 


Table 1. — Locations and Installation Dates of 
the Stations of the Apollo Seismic Network 


Station 

Date of 
Installation 

Location 

Latitude 

Longitude 

12 

Nov. 19, 1969 

3.04°S. 

23.42°W. 

14 

Feb. 5, 1971 

3.65°S. 

17.48°W. 

15 

July 31, 1971 

26.08°N. 

3.66°E. 

16 

April 21,1972 

8.97°S. 

15.51°E. 
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Owing to the extreme quiet of the lunar 
surface, the instruments can be operated at 
peak magnification of approximately 20 mil- 
lion. This is two to three orders of magnitude 
higher than can normally be used on Earth. 
All but two of the 16 separate seismometers 
are presently operating. The short-period 
seismometer at station 12 has failed to oper- 
ate since initial activation, and the long- 
period vertical component seismometer at 
station 14 became unstable after 1 yr of 
operation. 

The following section summarizes the ma- 
jor findings from earlier work for the bene- 
fit of the reader unfamiliar with these results. 
A recent review of the results has been 
presented by Lammlein et al. (ref. 1). 

Brief Review of Previous Results 

Prior to the first lunar landings, the seis- 
mic experiment team had anticipated that 
moonquakes, if extant, might originate along 
the large rill-like features that crisscross 
the lunar surface or along the fracture sys- 
tems associated with the large maria basins. 
High-frequency signals that are possibly 
shallow moonquakes have, in fact, been re- 
corded and are among the largest of the 
signals recorded to date. However, these 
high-frequency signals are extremely rare; 
11 such events have been identified over a pe- 
riod of 3 yr. Much more numerous, but very 
small in magnitude (maximum Richter mag- 
nitudes of about 2), are moonquakes that 
occur at great depth. These are concentrated 
in a relatively thin zone at depths of between 
about 600 and 800 km. A third category of 
moonquakes has been recognized. These are 
small events that can be detected at each sta- 
tion at maximum ranges of a few tens of 
kilometers. The frequency of occurrence of 
these events is closely correlated with the 
lunar surface temperature cycle. These are 
almost certainly of thermal origin. Slump- 
ing of material along steep slopes or disloca- 
tions along fracture or bedding planes of 
material very near the surface caused by 
thermal stresses have been suggested as pos- 


sible source mechanisms (ref. 2). These will 
not be discussed further here, except to re- 
mark that the mechanism leading to thermal 
moonquakes may represent an important ero- 
sional process on the lunar surface. The total 
annual energy release from moonquakes is 
estimated to be less than 10 15 ergs per year. 
This is about nine orders of magnitude lower 
than that of the Earth. If the sensitivities 
of the lunar seismographs were limited to the 
levels of their terrestrial counterparts, we 
would have recorded no moonquakes at all. 
Our conclusion would have been that the 
Moon is totally devoid of internal activity 
leading to quakes. Signals are recorded from 
several hundred meteoroid impacts per year 
at each station. By using the impacts from 
the third stages of the Saturn boosters as 
calibrations and a statistical method devel- 
oped by James Dorman (described in ref. 3), 
the present meteoroid flux in the vicinity of 
the Moon has been estimated to be 
log N = —1.62 — 1.16 log to 
where N is the cumulative number of me- 
teoroids of mass to (in grams) and greater 
which strike the Moon per year per square 
kilometer. This flux estimate is one to three 
orders of magnitude lower than those de- 
rived from earlier Earth-based measure- 
ments. However, it is in agreement with 
estimates derived from the distribution of 
crater sizes on the youngest lunar maria, 
assuming a plausible decrease in meteoroid 
flux between the time of maria formation and 
the present. Indeed, Wetherill (ref. 4) has 
pointed out that previous Earth-based flux 
estimates were too high by a factor of about 
35 to be consistent with the crater size dis- 
tributions observed on the maria. We esti- 
mate that a meteoroid of mass 7 to 10 kg 
can be detected by the Apollo Seismic Net- 
work from any point on the Moon. Thus, the 
Moon, when attached to a seismometer, has 
proven to be an unexpectedly good “sound- 
ing board” for the detection of meteoroid im- 
pacts. This results from a combination of 
two factors: the extremely high instrument 
sensitivities possible on the lunar surface 
and the unusually high efficiency of seismic 
wave propagation throughout most of the lu- 
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nar interior. The major weakness in the 
method thus far is the availability of only 
one calibration point (S-IVB) impact en- 
ergy). Obviously, a larger number of man- 
made impacts covering a wide range of 
impact energies would be highly desirable. 
Gault (personal communication) has pointed 
out that we may need to increase our flux es- 
timate by a factor of 2 to correct for the 
probability of nonnormal angles of impact. 
It should also be noted that attempts to use 
near-range seismic sources (thumper, etc.) 
as calibrations for meteoroid impact signals 
recorded by the SPZ component lead to 
higher meteoroid flux estimates (ref. 2). 
The discrepancy between these two methods 
remains to be resolved. 

Lunar seismograms differ markedly from 
typical Earth recordings. The most striking 
feature of lunar signals is their long dura- 
tion. Examples are shown in figure 1. Impact- 
generated lunar signals have emergent 
beginnings, increase gradually to a maxi- 
mum, and then slowly decay. Following the 
first one or two cycles of the P-wave, ground 
motion is very complex, with little or no cor- 
relation between any two components. The 
onset of the shear wave from near impacts 
(ranges less than 1000 km) cannot be iden- 
tified with certainty. Coherent surface wave 
trains displaying dispersion have not been 
recognized in any recordings to date, al- 
though scattered surface waves undoubtedly 
contribute to the signals. The direction of 
propagation can rarely be determined from 
the particle motion at a single station. Thus, 
many analysis techniques found useful in 
terrestrial seismology cannot be applied to 
lunar signals. Signals from the deep moon- 
quakes are generally less complex in the 
early part of the wave train than those from 
impacts. The largest amplitudes are clearly 
associated with the shear waves ; however, 
the exact onset of the shear waves is diffi- 
cult to identify. These unusual characteris- 
tics of lunar signals appear to be accounted 
for by assuming a high degree of hetero- 
geneity, resulting in intensive scattering, 
and very low absorption of seismic wave en- 
ergy in a thin surficial zone, referred to as 


the scattering zone. Nakamura (ref. 5) has 
successfully described the envelope charac- 
teristics of lunar signals by using diffusion 
theory. His analysis shows that most of the 
scattering occurs in the outer few hundred 
meters of the Moon in which the Q of the 
medium must range between about 3000 and 
5000. In fact, the extremely low absorption of 
seismic waves in this material is the domi- 
nant factor leading to the characteristics of 
lunar seismograms described above. Many 
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Figure 1 . — Compressed time-scale records of long- 
period (LP) components for two natural seismic 
events recorded at four seismic stations. X and Y 
refer to the horizontal component seismometers; 
Z refers to the vertical comvonent seismometers. 
Each successive trace deflection of an LP com- 
pressed time-scale record represents the sum of 
the absolute values of the differences between 40 
successive data samples during consecutive, non- 
overlapping 6-s time intervals. Successive deflec- 
tions are plotted on opposite sides of the zero base- 
line. The relative amplitudes of the signals are 
indicated by the brackets following each record. 
Tick marks are 10 min apart, and hours are la- 
beled. The category An moonquake zone is located 
at the eastern border of Mare Serenitatis. For 
category A moonquakes, the direct shear wave S is 
prominent on seismograms from the horizontal 
seismometers. The direct compressional wave P is 
observed on seismograms of the largest category A 
moonquakes. The category C meteoroid impacted 
roughly 500 km north of station 15. The onsets of 
the direct compressional and shear waves can 
rarely be identified on LP seismograms of category 
C events. 
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terrestrial seismograms would possess the 
same characteristics as those of lunar seis- 
mograms if the absorption of seismic energy 
were reduced to lunar values. There is a 
growing body of laboratory data supporting 
our earlier contention that this marked con- 
trast in absorption can be ascribed to the 
nearly complete absence of volatiles, princi- 
pally water, in the outer shell of the Moon. 

A seismic velocity model for the upper 150 
km of the Moon has previously been defined 
(refs. 3 and 6) on the basis of observed trav- 
eltimes and amplitudes of body waves from 
nine impacts of S-IVB and LM space ve- 
hicles, combined with data from laboratory 
measurements on returned lunar samples. 
Seismic signals from these impacts were re- 
corded at ranges from 9 km to 1750 km. 
Figure 2 shows a first-order model for the 
variation in the velocity of compressional 
waves, with depth consistent with the travel- 
times of seismic waves from these manmade 
impacts and a meteoroid impact that oc- 
curred within the array. The most important 
feature of this model is the abrupt increase 
in velocity in the depth range between 50 
and 55 km. By analogy with the Earth, we 
refer to the zone above this discontinuity as 
the crust and to the zone below as the man- 
tle. Near the surface the velocity increases 
rapidly from a value of about 100 m/s mea- 
sured in the regolith, reaching a value of be- 
tween 6.3 and 7.0 km/s at the base of the 
crust. The rapid increase in velocity near the 
surface can be explained by the progressive 
compaction of materials so thoroughly frag- 
mented that the velocity of propagation 
through them is determined primarily by its 
mechanical state rather than by its chemical 
composition. Velocities in the range from 6 
to 7 km/s are appropriate for, but not re- 
stricted to, the rock types found to predomi- 
nate in the lunar highlands (anorthositic 
gabbros and aluminous basalts) . At the man- 
tle boundary, the compressional wave veloc- 
ity increases to about 8.1 km/s — a value 
appropriate for rocks of olivine-pyroxene 
composition. As shown by Kovach and Wat- 
kins (ref. 7), the velocity variation in the 
outer crust may be stepwise in some areas 


where flows (ash or lava) have escaped com- 
plete obliteration, but the smooth variation 
shown in figure 2 appears to be a good ap- 
proximation of the stepwise increase they in- 
fer for the region of the Apollo 17 landing 
site. 



Figure 2 . — Compressional wave velocity model based 
on seismic data from nine manmade impacts and 
one near meteoroid impact. 
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It was found soon after the installation of 
station 12 that moonquake signals could be 
grouped into sets, with members of each set 
having waveforms that match one another in 
detail. Moonquakes in each group occur at 
regular intervals, normally once per month, 
and at specific times during the lunar orbit. 
In some cases as many as three or four moon- 
quakes of a given set will occur during a 
monthly cycle, but such multiple events oc- 
cur over a relatively short interval of a few 
days or less. The repetitious character of 
these events strongly supports the hypothe- 
sis that they are moonquakes, each group of 
matching signals corresponding to an active 
zone or focus within the Moon at which the 
repeating moonquakes originate. Forty-one 
sets of signals, corresponding to 41 active 
foci, have been identified to date. Since the 
total number of moonquakes in the matching 
groups represent only about 10 percent of 
the several thousand moonquakes recorded to 
date, it is quite probable that many other ac- 
tive foci exist. However, the signals recorded 
from them are too small to permit detailed 
analysis of their waveforms. 

Two tidal periodicities have been recog- 
nized in the moonquake activity patterns 
(variations in times-of-occurrence and am- 
plitudes): 13.6 days and 206 days. The 
shorter period variation corresponds to the 
nodical or draconic month. The longer period 
is introduced by solar perturbation of the 
lunar orbit, which results in a cyclical varia- 
tion of the Earth-Moon separation at times 
of perigee. A third periodicity, correspond- 
ing to the successive synchronizations of the 
anomalistic and nodical months and having 
a period of about 6 years, appears to be 
present. However, the sample length is not 
yet long enough to confirm this. It is evident 
that tides play a dominant role in the gener- 
ation of moonquakes. However, with very 
few possible exceptions, the polarities of seis- 
mic signals from a given source are identi- 
cal. This implies that the dislocation is 
progressive and not periodically reversed. 
Progressive dislocation suggests secular ac- 
cumulations of strain periodically released 
by moonquakes. Weak convection within the 


deep lunar interior or relaxation of tidal de- 
formation as the Moon recedes from the 
Earth have been suggested as possible 
sources of secular accumulation of strain. 

Using the velocity model described above, 
extrapolated to the depths of moonquakes, 
27 of the active foci have been located. The 
epicentral locations are shown in figure 3. In 
nine cases the depth of focus was assumed in 
the calculation. The moonquake foci are 
concentrated in two narrow belts. Both belts 
are 100 to 300 km wide, about 2000 km long, 
and 600 to 800 km deep. (Earlier estimates 
placed the depth of the moonquake zone at 
between 800 and 1000 km.) Whether the ap- 
parent alignment into belts is real or a mir- 
age that will disappear as more epicenters 
are located is not known. However, the dis- 
tribution is certainly not random. The 
epicenters lie approximately along arcs of 
great circles that intersect at an angle of 
about 80 degrees. This pattern cannot be ex- 
plained as a natural consequence of the dis- 
tribution of seismic stations, nor do we know 
of any tidal stress component that would 
lead to such a pattern. In addition, no sys- 
tematics in the relative times of occurrence 
of the moonquakes within a given belt have 
been discerned, except that peaks in the 
combined activity detected from all foci oc- 
cur at average intervals of 13.6 days and are 
approximately in phase at all stations. 

Only one moonquake epicenter has 
been located on the farside of the Moon. For 
signals from this source, it was noticed 
that shear waves were missing in the seis- 
mogram recorded at the most distant station 
for which the signal could be detected 
(station 14). Since shear waves from moon- 
quakes are normally prominent at station 14, 
this observation remained a mystery until 
the phenomenon of missing, or greatly de- 
layed, shear waves was observed in the 
records from a distant impact. Although 
other interpretations are possible, Nakamura 
et al. (ref. 8) found that these observations 
could be explained by assuming that partial 
melting begins in the present-day Moon at a 
depth of about 1000 km. (This boundary 
must be moved upward to a depth of about 
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Figure 3 . — Map of the nearside of the Moon locating the Apollo 12, H, 15, and 16 seismic stations (table 1) 
and the category A moonquake epicenters. Solid circles indicate the foci for which the depth can be de- 
termined. Open circles correspond to cases in which data are not sufficient for determination of depth. 
In these cases a depth of 800 km was assumed to locate the epicenter. Note that epicenters 1 and 6 and 
18 and 32 are so closely spaced that their separation cannot be distinguished at the scale plotted. One fo- 
cus occurs on the farside of the Moon. 
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800 km according to the more recent results 
presented here.) Following the same line of 
reasoning, significant melting cannot be 
present in the crust and mantle at depth shal- 
lower than the moonquake zone, since high- 
frequency shear waves are prominent in the 
seismograms from all nearside moonquakes. 
Since little or no seismic activity has been de- 
tected above the moonquake zone and sub- 
solidus temperatures are inferred, we refer 
to this dynamically inactive outer shell as 
the lithosphere and to the “weaker” central 
zone, in which partial melting is inferred, 
as the asthenosphere. 

Recent Results 

HIGH-FREQUENCY TELESEISMIC 
EVENTS 

A set of seismic signals with unusual char- 
acteristics has recently been identified (ref. 
9). They are distinguished from all other 
signals recorded at comparable ranges 
(greater than 1000 km) by their high- 
frequency content. Only 11 signals of this 
type, designated high-frequency teleseismic 
(HFT) signals, have been identified in the 
3-year period for which data are available. 
Some of these are among the largest signals 
detected. Recordings for one of these events 
are shown in figure 4, along with representa- 
tive signals from a moonquake and a meteroid 
impact. The relatively high-frequency con- 
tent of the HFT signal can be seen in this 
figure by comparing the amplitudes recorded 
by the long-period seismometers (LPX,Y,Z) 
with those recorded by the short-period seis- 
mometers (SPZ). The relatively large 
amplitude of the SPZ record for the HFT 
event indicates that most of the signal energy 
is at frequencies above 1 Hz; whereas fre- 
quencies below 1 Hz are dominant in signals 
for distant moonquakes and impacts. This 
comparison is made for a larger number of 
events in figure 5, where the ratio of the max- 
imum signal amplitudes recorded by the 
SPZ and LPY (one of the long-period hori- 
zontal components) seismometers is plotted 
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Figure 4 . — Seismogram of a high-frequency teleseis- 
mic (HFT) event detected on March IS, 1973 (cen- 
ter), compared with those of a moonquake detected 
on June 5, 1973 (left: category Ai), and a meteor- 
oid impact detected on May 13, 1972 (right). All 
of these seismograms were recorded at the Apollo 
16 station. LPX, LPY, and LPZ stand for three 
orthogonal components of a long-period instrument 
peaked at about 045 Hz. SPZ stands for a short- 
period vertical component peaked at 8 Hz. Note the 
relatively impulsive beginning of P- and S-wave 
arrivals of the HFT event, similar to those of the 
moonquake, and the large SPZ amplitude of the 
HFT event. 


IO r 


< 

oc 

LlI 

Q 

3 


CL 

5 

< 

>- 

Q- 


rsi 

Q_ 

(O 


1.0 


0.5 [ 


0.2 L 


o - HF EVENTS 
■ - MOONQUAKES 
a - METEOROID IMPACTS 
0 LM IMPACTS 


50 100 150 200 250 

P - S INTERVAL , s 


300 


Figure 5. — SPZ/LPY amplitude ratio versus S-P 
time interval for the HFT events and selected 
other events recorded by the Apollo H station. A 
large amplitude ratio means a greater high- 
frequency content of a signal. Increasing S-P time 
corresponds to increasing distance from the sta- 
tion. 
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versus the time interval between the P-wave 
and S-wave arrivals for each signal. This time 
interval is a function of distance between the 
source and the seismic station — increasing 
time interval corresponding to increasing dis- 
tance. The distinctive character of the HFT 
signals is evident in this plot. 

The beginnings of the P-wave and S-wave 
trains are much more impulsive for HFT 
signals than those of typical meteoroid im- 
pact signals, suggesting that scattering near 
the source is much less for the HFT events 
than for normal surface sources. On the other 
hand, the predominance of shear waves in the 
signals from the most distant HFT events 
requires that the source be shallower than 
300 km (ref. 9). 

Several source mechanisms for these 
unusual signals are suggested: either they 
are generated by meteoroids that penetrate 
the scattering zone into more competent 
material below (200 to 300 meters) , they are 
generated by impacts in areas where unusally 
competent and structurally homogeneous 
material lies close to the surface, or they are 
shallow moonquakes. The latter interpreta- 
tion is favored. It is difficult to conceive of 
any process that would preserve isolated 
patches of the lunar surface relatively un- 
disturbed by meteoroid impacts. If such 
regions do exist, gradational cases would be 
expected, resulting in a broad spectrum of 
signal characteristics instead of the sharp 
division between signal classes observed. On 
the other hand, the existence of some shal- 
low moonquake activity in a cooling planet 
would be expected. 

STRUCTURE AND STATE OF THE DEEP 
LUNAR INTERIOR 

Information on the structure and state of 
the deep lunar interior is derived principally 
by analysis of seismic signals from distant 
impacts and moonquakes. As noted above, 
meteoroid impact signals are so emergent 
that it is rarely possible to pick the onset of 
the P-wave with certainty, and the onset of 
shear waves can be separated out from the 


compressional wave train only in the signals 
from distant impacts. Signals from four of 
the largest impacts and two of the HFT 
events have been selected as suitable for 
the present analysis. The traveltime curves 
for these events, as given by Nakamura et al. 
(ref. 10), are shown in figure 6. These data 
can be inverted by the Wiechert-Herglotz 
method to give seismic velocities as a func- 
tion of depth, assuming that velocity gradi- 
ents are not so large as to invalidate the 
method. However, several points can be 
made by inspection of these curves without 
further analysis. First, we note that the 
curve for shear (S) waves departs markedly 
from the curve for compressional (P) waves 
at a range of about 85° (2600 km), corre- 
sponding to maximum depths for the S-wave 
ray path of about 300 to 400 km. The 
velocity of shear waves must begin to de- 
crease sharply (i.e., Poisson’s ratio begins to 
increase) in this depth range. Secondly, 
the P-wave traveltime for the most distant 
point is delayed by about 120 s from the value 
expected by smooth extrapolation of the 
curve. This indicates that the velocity of P- 



Figure 6. — Travel-time data for compressional (P) 
and shear (S) waves recorded from four distant 
impacts and two HFT events. 
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Figure 7. — Four-station expanded time scale play- 
outs of a category A, s moonquake detected at 22:53 
hours on January 5, 1973. X and Y refer to LP 
horizontal component seismometers; Z to the LP 
vertical component seismometers. Tick marks are 
1 min apart. These signals are typical of moon- 
quakes recorded at each of the stations. The direct 
shear wave (S) is prominent on the horizontal 
component seismograms at all of the stations. The 
direct compressional wave (P) is also observed at 
stations 14, 15, and 16. 


waves decreases abruptly at depths of be- 
tween 1380 and 1570 km. 

Below the zone of abrupt shear-wave 
decrease, the Wiechert-Herglotz method 
cannot be applied for shear waves with con- 
fidence, and signals from moonquakes that 
originate beneath this zone must be used. 
Signals from one of the moonquakes used in 
this analysis are shown in figure 7. Combining 
the traveltime data from deep moonquakes 
with the data of figure 6, the model shown 
schematically in figure 8 has been derived. 
For purposes of discussion, it is convenient 
to divide the lunar interior into five major 
zones: 


Zone 1 — The Crust 

Zone 1 represents a layer approximately 55 
km thick in the region of stations 12 and 14. 
If the crust is a global feature, as suggested 
by several lines of evidence (ref. 1), it is 
probably quite variable in thickness, with 
thickness of the backside crust substantially 
greater than that of the nearside (ref. 11). 


The crustal velocities are appropriate for the 
plagioclase-rich materials inferred from 
sample analysis, possibly grading from anor- 
thositic material at the top to more gabbroic 
material at depth. 


Zone 2 — The Upper Mantle 

Zone 2 is approximately 250 km thick, with 
a compressional wave velocity of about 8.1 
km/s at the top, decreasing to about 
7.8 km/s at the bottom. Possible composi- 
tions of this zone can be inferred from the 
velocities and densities corresponding to the 
model of figure 8. Based upon geochemical 
considerations (ref. 12), olivine and pyrox- 
ene are the most probable constituents of the 
rocks of the upper mantle. Curves of veloc- 
ity versus density for these minerals are 
given in figure 9, along with the values de- 
rived from our preliminary model, assuming 
a uniform composition for the upper mantle 
(i.e., the velocity decrease with depth in the 
upper mantle is due to increasing tempera- 
ture with depth). The experimental values 



, Figure 8 . — Pictorial representation ( drawn approxi- 
mately to scale) of the internal structure of the 
Moon as described in the text. X marks indicate 
the zone of deep moonquakes. 
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Figure 9 . — Density versus velocity (corrected to 
standard temperature and pressure) calculated 
from the seismic data of figure 6 for the upper 
mantle (zone 2) of the Moon. Compressional- and 
shear-wave data for the olivine (Fo — Fa) series 
(ref. 13) and the orthopyroxene (En — Fs) series 
(Chung, personal communication, 1974) are shown 
for comparison. The calculations permit variations 
in the various lattice dynamical parameters over 
ranges approximate for rocks of mafic composi- 
tion. Calculations have been made for two cases. 
(1) The crosshatched block corresponds to the as- 
sumption that the composition of the upper mantle 
is uniform throughout, i.e., that the effects of tem- 
perature and pressure balance to give the experi- 
mentally determined velocity-depth functions. This 
requires a velocity gradient of 4 to 5° C/km in the 
upper mantle, (2) The open block below shows the 
change introduced by decreasing the temperature 
gradient to between 2 and 3° C/km. For this case, 
the intrinsic velocity and density must decrease 
slightly with depth. From this comparison, a com- 
position close to the magnesian end of the olivine 
series (80 percent forsterite) is suggested for the 
upper mantle. For both cases described above, the 
effect of including an iron core ( zone 5) of density 
7.5 g/cc and a maximum radius of 360 km is also 
shown. 


for the upper mantle thus derived are in good 
agreement with the olivine data (80 to 85 
percent forsterite) . The calculated points are 
shifted slightly to the left if an iron-sulphide 
core (zone 5) is assumed. The assumption 
of constant composition requires a tempera- 
ture gradient of 4 to 5°C/km for the upper 
mantle. Alternatively, the traveltime data 
can be satisfied by assuming a variation 
from magnesium olivine at the top to more 
pyroxenitic material below and reduction of 
the temperature gradient to between 2 and 
3°C/km. Walker et al. (ref. 14) have 
pointed out that a mafic cumulate of thick- 
ness corresponding to that inferred here for 
the upper mantle (about 250 km) might be 
expected if a plagioclase-rich crust 50 to 60 
km thick is present. 

Zone 3 — The Middle Mantle 

This zone, extending from about 300 to 
800 km, is characterized chiefly by the re- 
duced velocity of shear waves (3.6 to 4.0 
km/s) within it. It is likely that the com- 
pressional wave velocity is also decreased 
slightly in this layer, but this cannot be re- 
solved with the available data. Thus, the 
material of zone 3 must have a high Poisson’s 
ratio (0.34 to 0.35) compared to that of the 
material above in which the Poisson’s ratio 
is approximately 0.25. One possible inter- 
pretation is that this zone represents primi- 
tive lunar material below the zone of initial 
melting that produced zones 1 and 2. Mare 
basalts may have been produced by partial 
melting within zone 3, 3.2 to 3.8 b.y. ago in 
accordance with the proposals of Duba and 
Ringwood (ref. 12) and others. Moonquakes 
are concentrated within a relatively thin 
zone at the base of zone 3, in a depth range 
of from 600 to 800 km. Because high- 
frequency shear waves are recorded from all 
of these sources (except the farside focus), 
we can say immediately that widespread 
melting in zones 1, 2, and 3 is not possible, 
i.e., temperatures in these zones must be sub- 
solidus. In addition, there is little or no scat- 
tering of seismic waves in zone 3. Thus, if 
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this represents an accumulation of primitive 
materials, either individual blocks are small 
compared to a wavelength (less than about 
1 km in dimension) , the elastic properties of 
the individual blocks do not vary appreciably, 
or any inhomogeneities have been removed 
by subsequent melting. 

There is some evidence that the transition 
from zones 2 to 3 may be relatively sharp. 
A prominent phase, designated R in figure 
10, has been identified in the wavetrains of 
the largest moonquake signals. The relative 
arrival times for this phase can be explained 
by assuming that the phase is generated by 
conversion from S to P at an interface 
located at a depth of about 300 km, the 
approximate depth for the boundary between 
zones 2 and 3 derived from shear-wave travel- 
time data. Similar mode conversions have 
been tentatively identified in records from 
deep earthquakes (ref. 16). 

Zone A — The Lower Mantle 

This zone is characterized chiefly by the 
absence of identifiable shear waves for events 
in the distance range for which shear waves 
would enter this zone. This can be ascribed 
to (1) a sharp decrease in shear-wave 
velocity at the top of the zone, leading to a 
shadow zone for shear waves, or (2) high 
attentuation of shear waves within this zone. 


The latter interpretation is favored for the 
following reasons. First, we must explain the 
concentration of moonquake activity, closely 
correlated with lunar tides, immediately 
above this zone. Calculation of tidal energy 
density as a function of depth shows that a 
sharp maximum is introduced in the depth 
range of moonquakes if a zone of much 
“weaker” material (zone 4) is introduced be- 
low (ref. 16). Second, the velocity of P- 
waves cannot decrease more than 0.2 to 0.3 
km/s within this zone. A greater decrease in 
compressional wave velocity would be ex- 
pected to accompany the sharp decrease in 
the shear-wave velocity required to produce 
an optical shadow zone. To explain these ob- 
servations, Nakamura et al. (ref. 8) have 
argued that partial melting is probable in 
the lunar mantle. If so, a temperature of ap- 
proximately 1500°C is inferred for the top of 
zone 4, assuming a mafic composition. 

Zone 5 — The “Core” 

This zone is characterized by a sharp de- 
crease in the velocity of P-waves beginning 
at a depth of from 1380 to 1570 km. The 
inferred P-wave velocity for this zone is 
about 5 km/s, suggesting that the material 
of this zone is completely molten. Unless the 
material of zone 5 is more dense than that of 
zone 4, convective instability between zones 
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Figure 10 . — Seismograms from the long-period vertical component for two moonquakes, showing a phase 
(R) that may correspond to conversion from shear wave (S) energy to compressional wave (P) energy at 
an interface at a depth of about 300 km (the boundary between zones 2 and S). 
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4 and 5 would result. It is suggested that 
zone 5 represents a molten core of iron sul- 
phide in accordance with the models earlier 
suggested on geochemical grounds by Brett 
(ref. 17) and Duba and Ringwood (ref. 12). 
If the composition of the core is iron sulphide, 
it represents less than 1 percent of the total 
mass of the Moon and falls easily within 
moment-of-inertia constraints. However, this 
interpretation must be regarded as quite 
tentative pending the acquisition of more 
seismic data pertinent to the properties of 
this zone. 

Conclusions 

An iron sulphide core of radius 200 to 300 
km, as suggested by present seismic data, 
may have been the source of an early dipole 
field that magnetized lunar surface rocks. 

If accretional energy were responsible for 
melting the outer shell of the Moon to a 
depth of about 300 to 350 km to produce the 
crust and upper mantle, very short accretion 
times are required — probably less than 
10 000 years (ref. 12). 

The low level of lunar seismicity and its 
concentration at great depth compared to the 
Earth now appear to be explained by the 
differences in thickness of the lithospheric 
shells and internal thermal energy of these 
two planets. The lunar lithosphere is simply 
too thick (about 800 km) to be disrupted by 
the relatively weak convective motion that 
might be induced by temperature gradients 
below. In the absence of this source of seismic 
activity, i.e., collisions between lithospheric 
plates, dislocations induced by tidal defor- 
mation of the Moon, are the dominant mech- 
anisms for generation of quakes. Assuming 
that the Moon is cooling (contracting) 
slightly at present, as suggested by the 
thermal models of Toksoz and Solomon (ref. 
18) , the accumulation of compressional stress 
in the lithosphere would be expected to result 
in moonquakes at shallow depth. The few, 
relatively energetic events of unusually high- 
frequency (HFT) signals that have been 
detected may be explained in this way. 
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